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We analyze the systematic errors in the positions of bursts in the 
BATSE IB, 2B and 3B catalogs, using a likelihood approach. We 
use the BATSE data in conjunction with 196 single IPN arcs. We 
assume circular Gaussian errors, and that the total error is the sum in 
quadrature of the systematic error <7 sys and statistical error cr sta t, as 
prescribed by the BATSE catalog. We find that the 3B burst positions 
are inconsistent with the value <r sys = 1.6° stated in the BATSE 3B 
catalog. 
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INTRODUCTION 



■ The stated systematic error in the BATSE 3B catalog (|l|) burst locations 
, is (T S ys = 1.6°. This value was estimated by taking the RMS deviation of 

BATSE positions from known positions of 36 bursts, determined by the IPN, 
WATCH, and COMPTEL (g). Unfortunately, many of these same known 

■ positions were used to calibrate the BATSE burst location software, i.e., as a 
I , guide in determining what effects to include in the burst location algorithm. In 

i-H ' order to calibrate properly the BATSE burst position errors, an independent 
<0> ', set of burst locations is needed. Such a set exists, in the form of 196 bursts 
for which single IPN arcs exist (||). In this paper we analyze the systematic 
\ errors in the positions of bursts in the BATSE IB, 2B and 3B catalogs using 

■ these 196 bursts. Our analysis is based upon the likelihood approach. 
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CZ2 . We assume that the systematic error a sys and statistical error er s t a t in po 
^ ' sition are circular Gaussians, to be added in quadrature, as prescribed by the 
BATSE catalog. The circular approximation should be good, since roughly 
2/3 — 3/4 of all BATSE bursts have \ 2 positional contours that are nearly cir- 
cular (||), and the bursts with IPN arcs tend to be those with larger fluences. 
C\3 We further assume that the sky is flat (i.e., that the size of the total error 
in the BATSE burst position is not too large), and that the IPN arcs have 
zero width. The former is a good approximation, since the largest BATSE 
a tot ~ 12° <C 90° among the bursts with IPN arcs. The latter is a good ap- 
proximation since the characteristic width of the 3-er contours of the IPN arcs 
is a few arcminutes, which is much less than the smallest BATSE 0t o t ~ 1°. 
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Figure 1. Configurations for a burst with a known position (left diagram) 
and for a burst with a single IPN arc (right diagram). 

Figure 2. Probability density as a function of a sys for the ML constant cr sys 
model. The vertical lines show the ±l-cr interval. 

The likelihood function for bursts with known positions (e.g., two intersect- 
ing IPN arcs) is given by 

N N i ^ d 2 

£ (CTsys )=n^-riw exp (-2#' (i) 

i— 1 i— 1 ^ ** 

whereas the likelihood function for bursts with single IPN arcs is given by 

JVarc N arc 2 

Here, rfj is the deviation of the BATSE position from the known position, Oj 
is the perpendicular deviation of the BATSE position from the IPN arc (see 
Figure 1), and m = 0.43a to t = 0.43(a 2 ys + (T 2 tat ) 1/2 . 

The likelihood function C(a sys ) allows exploration of any model that we 
may have for a syB . Here we focus on three models: (1) cr sys = a constant, 
(2) a sys = A{S/10- r ° erg cm" 2 )", and (3) <7 sys = ^(cr st at/l°) Q . The first 
model has one free parameter (<7 sys ); the second and third models have two 
free parameters (scale factor A and power law index a). 



RESULTS 

Fitting the constant cr sys model to the 18 BATSE 3B bursts with IPN 
positions, we find a maximum likelihood (ML) value cr^S" = 1.85°1;qooo. This 
is consistent with the value of 1.6° quoted in the BATSE 3B catalog (Q), which 
was found using these bursts and some others. 

Fitting the constant er sys model to the 196 BATSE 3B bursts with IPN arcs, 
we find a ML value o^J = 3.7°to?,f°- Tnis is inconsistent (at the 10er level!) 
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Normalized Deviation Normalized Deviation 

Figure 3. Comparison of expected and observed cumulative distributions of 
the normalized deviation D = \ai/a\ ot \ for the burst samples and data as 
labeled, assuming the ML models of a sys as labeled. 

with the value of 1.6° quoted in the BATSE 3B catalog, as shown by the 
probability distribution in a sys (sec Figure 2) and comparison of the expected 
and observed distribution of normalized deviations D = \ai/a\ ot \ assuming 
<7 sys = 1.6° (see Figure 3). 
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Fitting the second model, in which tr sys = A(S/10~ 5 erg cm -2 )", we find 
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Figure 4. Upper left-hand and lower and upper right-hand panels: 2-cr, 
and 3-er contours in the (A, a)-plane for the burst samples and data as labeled, 
assuming the cr st at-dependent model. Lower left-hand panel: same as Figure 
2, except for the cr s t a t-depcndent model. 



ML values A ML = 3.8 ± 0.3 and a M L = -0.18+^|. The probability P(a > 
0) = 1.6 x 10 -3 , showing that the data prefer the power-law model over the 
ML model with constant (T sys = 3.7°. 

Fitting the third model, in which a sys — Aa" tat , we find ML values Aml = 
3.8 ± 0.2 and a ML = 0.321°;°^ (see Figure 4). The probability P(a < 0) = 
4.7 x 10 , showing that the data strongly prefer the power-law model over 
the ML model with constant er sys = 3.7°. 

These results imply that a sys is strongly correlated with S and <7 s t a t- For 
example, cr S ys(CTstat = 0.1°) = 2° while cr sys (cr s tat = 10°) = 8°. This is il- 
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lustrated in Figure 5, which shows the ML S-dependent and cr sta t-dcpcndcnt 
models. 
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Figure 5. Left-hand panel: The ML S-dependent model. Right-hand panel: 
The ML <T s tat-dependcnt model. The points are not data points, but rather 
are the ML constant a sys model for the given interval in S or cr s t a t ■ 

It is also interesting to examine the character of (T sys for subsets of the 3B 
catalog, in particular the IB sample of bursts in which evidence for repeating 
was found. Fitting the three models to the "old" IB positions, we find a sys — 
4.0° for the ML constant <j sys model, and that the data do not request a more 
complicated model(see Figures 3 and 4). In contrast, fitting the 3B positions 
for the IB sample of bursts, we find that the data request the cr s t a t-depcndent 
model at the 95% confidence level [e.g., P(a < 0) = 1.6 x 1CT 2 ] (again sec 
Figures 3 and 4). 

These results imply that for the IB sample of bursts, the new (3B) cr sys 's 
are smaller than the old (IB) ones for bursts with S > 4 x 10 -6 erg cm~ 2 or 
with Ustat < 1.6°. However, for the s=s 80% of bursts with lower fluences or 
larger <7 s t a t's, they imply that the new <7 sys 's are larger than the old ones. 

We are grateful to Kevin Hurley for providing us with the 3rd catalog of 
IPN arcs, without which this analysis would not have been possible. We also 
wish to thank Geoff Pendleton and Michael Briggs for discussing BATSE burst 
location issues with us. 

REFERENCES 

1. C. A. Meegan, et al., electronic catalog (grossc.gsfc.nasa.gov). 

2. G. N. Pendleton, private communication (1995). 

3. K. Hurley, et al, Third IPN Catalog (1995). 



